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I. INTRODUCTION 



By now, the Casimir force between a sphere and a plane was calculated using the functional determinant 
(or TGTG- or T-matrix-) representation and, for scalar field, numerically by the world line methods, 
and, of course, in the Proximity Force Approximation (PFA) . The extension of these methods to finite 
temperature is in its very beginning. Even for the PFA we did not find in literature a representation 
at finite temperature, although the calculation is quite simple. Recently (T], using the world line 
method, the temperature dependence was studied and an interesting interplay between geometry and 
temperature was observed. Using the functional determinant representation in j| a partly numerical, 
partly analytical study at medium and large separation was represented where also dielectric materials 
were included. It must be mentioned that the corrections beyond PFA are a topic of actual interest 
also for high precision measurements of the Casimir force where it is necessary to know the temperature 
corrections at small separation. 

In this paper we calculate the free energy at finite temperature for a sphere in front of a plane for 
conductor boundary conditions using the functional determinant method which gives the exact result. 
We pay special attention to the case of small separation and to the limiting cases of low and high 
temperature. Also, for small separation we establish the relation to the Proximity Force Approximation 
(PFA) and determine where it is applicable. 

The paper is organized as follows. In the next section we prepare the necessary formulas of the 
functional determinant representation. In sections 3 and 4 we calculate the temperature dependent part 
of the free energy at small separation by both methods, PFA and the functional determinant method. 
In sections 5 and 6 we derive the low and high temperature expansions using the functional determinant 
method. Finally, the results are discussed in section 7. 
Throughout the paper we use units with H = c = &b = 1. 



II. BASIC FORMULAS 



In this section we collect the formulas we need for the free energy in sphere-plane geometry. For the 
free energy we use the standard Matsubara-representation. For the interaction energy between a sphere 
and a plane we use it together with the functional determinant representation, 

T= \ E Trln(l-M(U), (1) 

n— — oo 

where £„ = 2nTn are the Matsubara frequencies. At zero temperature, the Matsubara sum becomes an 
integration, T X)^L_oo ~~ * 2 with £„ — > £ and the free energy turns into the vacuum energy, 

1 f°° d£ 

Eo= 9 ^Trln(l-M(£)). (2) 

In these formulas, M(£) is a matrix in the orbital momentum indices I and V , 

I l+V 

M l)l ,(o = d l mJ^ z E Kv»m)Hi; . (3) 
* ? i"=\i-i'\ 

For Dirichlet boundary conditions on the sphere, 

results from the scattering T-matrix. In these formulas, // +1 / 2 (:e) and K l+1 / 2 (x) are the modified Bessel 
functions. The geometry is shown in Fig. 1. We introduced the notations v = I + 1/2, v' = I' + 1/2 and 
v" = I" + 1/2, which will be used throughout the paper. 

The above formulas can be found in the original papers on the functional determinant method [Ej]. 
Here we follow the notations used in Q and in Q. The factors H\ v in ([3]) result from the translation 
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formulas. Their explicit form is 

H%, = v/(2Z + l)(2Z' + l)(2/" + 1) 
I V I" \ { I V I" 



m -m 



(5) 



where the parentheses denote the 3j-symbols. For Neumann boundary conditions on the sphere we have 
to substitute di(x) by 

N _ (MflM)' 

(svio/v?) 

and for Neumann boundary conditions on the plane we have to reverse the sign in the logarithm in ([1]). 

For the electromagnetic field the matrix M acquires indices for the two polarizations, which correspond 
to the TE and the TM modes in spherical geometry, 

, i+i' 

Ml ' v = \ik £ Kv»(XL)H l u, (7) 
» ^ i"=\i-y\ 

/ai;;, hAfdrm o 
1 a m , Af;, J 1, o -d™m 



with the factors 



v/Z(Z + l)Z'(/' + l) 



A' 



A«/ = , 2 "* L (8) 

which follow from the translation formulas for the vector field. The factors resulting from the scattering 
T-matrices are 

which is the same as for the scalar field with Dirichlet boundary conditions, and, for the TM-mode, 

dT u m = iumvm^ (10) 

When inserting these expressions into ([TJ or ((2J, the trace must be taken also over the polarizations. 
More details are given in Q . 



4 



III. THE PFA AT FINITE TEMPERATURE 



In order to discuss the free energy at small separation which is the relevant scale for PFA we first 
define the corresponding temperature regions of interest. These are 

1. low temperature: dT < RT < 1 , 

2. medium temperature: dT < 1 < RT , (11) 

3. high temperature: 1 < dT < RT . 

In each case d <C R holds. Further we introduce 

~3 (12) 

as small parameter. 

Below we consider the PFA at finite temperature in some more detail as in section 5.1.2. We do this 
for the electromagnetic field with conductor boundary conditions. For that, first we need to remember 
the free energy .Fn per unit area for two parallel planes. It can be obtained in several ways, for instance 
as special case from the general formula (JXJ) 

T ll (d,T) = T± J *L l n (WW^), (13) 



n— — 00 



where k is the momentum in parallel to the planes and are the Matsubara frequencies. By representing 
this free energy in the form 

F\\(d,T) = -^(l + g(2dT)) (14) 

(2d = 1/Tcff is the inverse effective temperature) we defined the function g(x) which describes the 
temperature dependence. From (|13p it can be written in the form 

g(a;) = _i + !^^_^g f°dfcfcin(i- e -V(2— rwy (15) 

A sum representation can be obtained from (|13[) by expanding the logarithm into a series with subsequent 
integration over k and summation over n, 



g(x) = -1 +45x 4 



m— 1 



coth(m7ra) . . ,._2 

— h (rmrx smh(m7ra;)) 



(16) 



Another representation can be obtained from (fT5|) using the Abel-Plana formula, see eq. (7.81) in [7] for 
example. Comparing that representation with (|16[) the inversion symmetry 

g(x)=x±g(±), (17) 

which holds for the function g(x), can be established. The asymptotic expansions of this function for 
small and large arguments are 



45C(3) 



X — 1 + . . . for x — > 00. 



The upper line follows directly from (fTo]) and the lower line from (IT71) . In (TT5)l . the dots denote expo- 
nentially small contributions. The asymptotic formulas reproduce g(x) quite accurate everywhere. Even 
at the symmetry point x — 1 the relative deviation is less than 5%. 

Now we apply the idea of PFA to (fT4"l) . We have to integrate T\ \ (d, T) over the area of a circle of radius 
R in the (x, y)-plane, 

J" PFA = I ' dxdyF\\{d + h(x,y),T), (19) 
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where d + h(x, y) is the separation between the plane and the sphere at the point (x, y) in the plane. 
Using polar coordinates we get 



J" PFA = / drr I d<pF\\{d + R{l- cos 0),T) 
Jo Jo 



(20) 



with r = Rs'mO. Using the azimuthal symmetry and changing for the variable t = 1 — cos 9 we come to 
the representation 



./o 



The corresponding approximation for the force, 



(21) 



rPFA _ x-PFA 

1 3d 



can be written in the form 



/ PFA = 2ttR (d,T) -J dt^d + Rt, 



T) 



(22) 



(23) 



In deriving (|2"3"|) . for the derivative of T\\ we used d/dd = {X/R)d/dt and integrated by parts. 

Both expressions, (|2T|) and (|23|) . are meaningful for e — > only. However, further simplification 
depends on the specific properties of the functions involved. 

The simplest case is T — 0. Here the function 0) is just the vacuum energy 



Eo 



720d 3 ' 



It increases proportional to 1 / <i 3 with decreasing separation d. In the limit e 
the parenthesis in (|23[) does not contribute and one arrives at the simple rule 



/ | P T F ^ = 2^^||(d,0) 



(24) 

0, the second term in 
(25) 



called proximity force theorem in [Iol | . 

For a general function F\\(d,T), the validity of (f25j) depends on the behavior T\\{d, T) ~ for 
d — > 0. From Eg. (1231) it is easy to show the following. If a > 0, i.e., if J-\\{d,T) increases for decreasing 
separation, the integral term in (|23p is by a factor d (or ci In d for a = 1) smaller than the first term in 
the parenthesis in (|2"3"|) and (|23|) holds. Contrary, for a < the integral term is dominating and (|2~5j) 
does not hold. A constant contribution to T\\{d,T), having a — 0, does not contribute to the force and 
([25]) does not hold either. 

Now we consider (12H and (|23|) with the free energy pT|) of parallel planes inserted. We get 



.F PFA = - 



720d 2 



1 



1 + e 



2e^ 



1 

1 j (e + *)3 



for the free energy and 

/P FA 



TT 3 R 

360d 3 



1 + 3/2e 
(1+e) 2 



,g(2dT) 



1 dt g(2(e + t)RT) 



(s + tf 



(26) 



(27) 



for the force. These expressions are meaningful in the limit of e — > only. First of all we mention that 
the corrections to the vacuum energy in the first terms in the square brackets in both expressions, which 
are due to the geometry, are known to be beyond the precision of PFA and we drop them in the following. 
In order to perform this limit in the temperature dependent parts it is meaningful to consider two cases. 
First we consider low and medium temperature as defined in (fTTj) . For this we substitute 



dT = sRT 



(28) 



G 



in (|26p and (j27|) and consider the limit s —> with i?T fixed. In the integrals we can put e = directly 
and using (fT5)l we obtain 



^-PFA = _ 



7T 3 i? 

720d 2 



1 + 2e 2 [ dt (1 - t) 
Jo 



g{2tRT) 
t 3 



for the free energy and 

/P FA 



TT 3 R 

' 360d 3 



l + £ 3^2 ( *T)3- £ ; 



1 _ u g{2tRT) 
a t 3 



(29) 



(30) 



for the force. These expansions are uniform in RT. 

We note that the temperature correction to the free energy does not depend on the separation d and 
that only the next order contribution not shown in (|29p gives the contribution to the force which is shown 
in (f30|). In opposite to (fT8|) , in the above formulas the dots denote contributions which are suppressed 
only by powers of the small parameter. 

As special cases, from (|2T)|) and we obtain for dT <C 1, i.e., for low temperature, 



T 



PFA 



■K A R 

720d 2 



1 + £ ^(*T)3 + 



for the free energy and 



PFA 



TT 3 R 



[l + e 3 8(i?T) 4 



(31) 



(32) 



360d 3 

for the force. The other special case is RT ^> 1. It corresponds to medium temperature and we come to 

n 3 R 



for the free energy and 



PFA 



PFA 



720d 2 



ir 3 R 
360d 3 



[1 + e 2 20(RT) 2 



r3 360C(3) (i?r) 3 



(33) 



(34) 



for the force. In (|33|) we used that the integral J °° dtt 3 g(t) = 5/2 could be calculated explicitly. 
Now we consider medium and high temperature as defined in (1111) . For this we substitute 



RT = 



dT 



(35) 



in (|26|) and (f27| and consider the limit e — > with dT fixed. In this limit, (|26|) and (|27|) turn into 

7T 3 i? 



Jr PFA = 



720d 2 



[1 + fc(dT) 



(36) 



for the free energy and 



/PFA = _ 



360d 3 



[1 + g{2dT) 



for the force. In (l36l) we defined 



h(x) = 2 



dt 



The limiting values of this function are 

h(x) ~ | 



20 x 2 + . 



g(2to) 
t 3 ' 

for x — > 0, 
for x — ^ oo, 



(37) 



(38) 



(39) 
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where the dots denote exponentially small contributions. We will meet this function again in section 4. 
We mention that (|22p holds for the contributions displayed in (|3"r?)) and (f3"T|). see Eq. ([70| . 

The expansions (l36l) and (|37|) are uniform in cff. For dT <C 1 these turn just into (|33|) and (|34|) , 
correspondingly. In this way, at medium temperature, the asymptotics of (|33j) and (|34[) for RT — > oo 
just match the asymptotics of (|36j) and (|37|) for dT — > 0. For dT ^> 1, i.e., for high temperature, useing 
the lower line in (I3TH) the asymptotics are 



for the free energy and 



J .p FA= _c(3p: + (40) 



PFA ((3)RT 
J 4d 2 ' [ ' 



for the force. Here the dots denote exponentially fast decreasing contributions. 

In this way, with equations ([!?§)) and (I3U1) we derived the free energy and the force in PFA for small and 
medium temperature. For fixed RT these are small, suppressed by 2 resp. 3 orders of e. Obviously, the 
rule (|25|) docs not hold. The square bracket in (|30|) should be equal to the parenthesis in (|T4|) taken for 
small dT but as seen from (TT5|) it is different. This is in agreement with the general discussion following 
the rule (|23|) since for small dT the thermal contribution to J-\\ has a < 0. It is only in the limiting case, 
for RT > 1, Eq.®, that (J25J) holds. 

For medium and large temperature the free energy and the force are given by Eqs. (|36[) and (I37|) . In this 
case the thermal corrections are not suppressed by powers of s. Their seize depends on dT. For small dT 
they are small and match (|29p and ([3"U|) taken at large RT. For dT ~ 1 they are of the same order as the 
vacuum energy and for dT ^> 1 they dominate. For the force (|3T|) . the rule (|2"5j) holds since in this case 
one has to take the function g{x) for large argument in the parenthesis in (fT4"|). This is also in agreement 
with the general discussion following the rule (f25j) since this behavior corresponds to an a > 0. In this 
way, the rule (1251) can be generalized to finite T in the region of medium and large temperature but not 
at small temperature. It must be mentioned that at small temperature the temperature dependent part 
is very small and therefore not relevant in applications. 



IV. THE FREE ENERGY AT SMALL SEPARATION 



In this section we consider the exact expression (TTJ for the free energy at small separation, i.e., for 
e <C 1. First we assume fixed values of R and T. In the sense of the definition (TTT1) this is the region 
of low and medium temperature. Also we restrict ourselves to the case of a scalar field with Dirichlet 
boundary conditions. We make use of a number of notations and the methods developed in [1, [111. We 
start with expanding the logarithm in (fTJ, i.e., with the representation 

?=-\ £ Ejtt £ £ II £ H 60 ' (42) 

n=— oo s=0 ro=— oo;=|m| Vj=l n 3 =|m|— I J 

where we defined 

s 

Z(0 =nM«+„ 1 ,, +Tl4+l (0 (43) 

i=0 

for the product of the (s+1) matrices M(£). We adopt the formal setting no = n s +i = 0. Up to 
the frequency sum, this representation is identical to the corresponding one in [6j. For instance, the 
convergence properties are the same. For any finite e, all sums in (|42l) converge. However, for e becoming 
smaller, a higher and higher number of terms give significant contributions until in the limit e — > the 
convergence gets lost. 

In [6|, the behavior for e —> was established by making an asymptotic expansion. First, the orbital 
momentum sums were substituted by corresponding integrations, afterwards a substitution of variables 
allowed for performing the limit. We will use that later in this section but first we proceed in a slightly 
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different way. We first transform the Matsubara sum into integrations using the Abei-Piana formula. In 
this way we come from ((32) to 

T = E + T T (44) 

with 

^ = -5EtJt/ i^E E II E *(*(*)-*(-*))• (45) 

s=0 J -°° m=-ool = \ m \ \j=l n s = \m\-l J 

The first term in the rhs. of (|44j) is just the vacuum energy ([2]) which appeared from the direct substitution 
of the Matsubara sum by the integration over £. The second contribution, Tt, involves the Boltzmann 
factor, 

^)= e | € |/r_ r (46) 

and it can be viewed as the 'pure temperature contribution'. The analytic continuation £ — > i£ is, so 
to say, the return to real frequencies. From the structure of the integrand one can show that there are 
no singularities the integration path could cross during rotation. Especially there are no singularities on 
the real frequency axis. If there were some these would correspond to discrete eigenvalues of the Laplace 
operator in the region between the sphere and the plane. However, since that is an open geometry, the 
spectrum is pure continuous. 

Now, for the vacuum energy Eq we can use the known asymptotic expansion for e — > 0. In leading 
order it is the PFA. Including the first correction beyond PFA it was calculated in [6] for the scalar field 
and in Q for the electromagnetic field. For example, for Dirichlet boundary conditions on the sphere 
and on the plane it reads 

C(4) 1 A , 1 



E » = - ^{ 1+ r + -)- (47) 

For the other boundary conditions this formula holds too, only the coefficients change. For the electro- 
magnetic field logarithmic contributions appear. 

Now we consider the pure temperature contribution (|45p . In order to perform the analytic continuation 
£ — » ±i£, we use the corres pon ding analytic continuations of the Bessel function entering ([3]) and (|4|). 
These are (see, for example fl2j ]) 

I v {iz) = i"J v {z), K v (iz) = ^C v H ( ?\z), (48) 

I v {-iz) = r v J v {z), K v {-iz) = % ^i v Hl 1 \z). (49) 

We insert these into the matrix M(£), <j3j> , and obtain 

Mi ^ ) = 4^\fm £ (-i) ( ' + ''-'" )/2 ^( 2 ei)^:, ( 50) 

where we took into account that only I" contribute for which (I + I' — I") is an even number. Next we 

( 1 2) 

use Hi, ' (z) — J v \z) ± iY v {z) in order to separate the imaginary part, 



r l + r u (ZR)r v „(2£L) 

2 



. . rv(eR)-rv,(2£L) 

1 — 2- 



l + r v (£R)rv»(2€L) 
where we defined the ratio 



(51) 
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TABLE I. The temperature dependent part of the force from the exact formula and from the PFA for e = 0.01 
and T = 1. 



R 


0.5 


1 


3 


h 


0.14 


0.59 


5.1 


/•PPA 
JT 


0.047 


0.33 


4.6 



In fact we need the imaginary part of Z{i£), i.e., of the product of (s+1) matrices (|5"Tj) . We observe, 
that the imaginary part of this product always contains at least one factor r„ (£) . 

Now we consider the convergence properties of the temperature dependent part (|45|) . In opposite to 
the zero temperature case, the convergence of the integration over £ is now provided by the Boltzmann 
factor. We are left with the question about the convergence of the orbital momentum sums. It is known, 
that these cease to converge for e — > in the vacuum energy where the modified Bcsscl functions enter. In 
the temperature dependent part we have the 'non-modified' Bessel functions. We need their asymptotic 
expansion for large indexes. Again using [l2j we note 

J v{z) = yilI e -Wn„ + l + „l„(,/2) (^ + (A\ Yu {z) = e W^-l-,ln(,/2) ^ + Q fA^ 

From here it follows that the ratio rv(£) is fast decreasing for increasing index v. As a consequence, 
the imaginary part of Z entering Tt, makes the orbital momentum sums convergent irrespective of e. 
We note that this is in opposite to the real part where the convergence comes from the combination 
J v i (£R)Y v >t (2£L)/Y V (£R) and is present only if R < L, or equivalently, e > 0, hold. 

In this way, the integral and the sums in (|4"51 remain finite if we put e — directly in Z. In this 
way, J~t has a finite limit for e — > 0, RT fixed. We mention that the same result was obtained in the 
preceding section using PFA. Therefore, from the representation of the free energy used in this section, 
which is an exact one, we conclude that the temperature dependent part of the free energy is at small 
separation a next-to-next-to-leading order correction. The leading order is given by the vacuum energy. 
It is proportional to 1/e 2 . After that there come the beyond leading order corrections proportional to 
1/e to the vacuum energy which were calculated analytically in 0, § [ll|. Only after that, in order e°, 
the temperature corrections start to contribute. 

The above discussion holds for R and T kept fixed and d small. This covers the region of low temper- 
ature until the border between low and medium temperature with RT ~ 1. In this region a numerical 
evaluation of the temperature dependent part is possible. We have made this calculation using the 
representation 

•Ft — ~z / ^ni(OiTr[ln(l-M(^T))-ln(l-MHeT'))] ) (54) 
27T Jo 

which is the same as Eg. (1451) but without expanding the logarithm and after the substitution £ — > £T. 
The results are represented in Fig. [5] and [3] for T = 1 as function of R. Fig. [5] is for e — and Fig. [3] for 
e = 0.1. For small R the calculation is quite easy since only the lowest orbital momenta contribute. In 
the limit of R — > the temperature part Tt is the same as obtained from the low temperature expansion 
which will be considered in the next section. For increasing R, more and more orbital momenta need 
to be included. In the process of evaluation we increased this number until the relative change dropped 
below 10 -2 . For example, for R = 2 we had to include orbital momenta until I — 9. 

An interesting feature of this calculation is that it allows for a direct comparison with the PFA results 
derived in the preceding section. The corresponding values for J 7 ^ calculated from (1211) are shown in 
the Fig. [5] and UJ as dashed lines. It can be clearly seen that the PFA deviates from the true values 
of Tt quite significantly. This statement holds for small temperature and at least up to the border 
between low and medium temperatures. For example, in Fig. [3] we cover the region from R = until 
R ~ 6 which by means of the definition (jlip is the beginning of the region of medium temperature. 
At R = 6, the PFA result for Tt is only about 70% of the true value. The picture becomes better for 
PFA if calculating the temperature dependent contribution to the force by means of (|22|) . We made a 
calculation for several values of e and R with T = 1 . The results are displayed in Tables U and [TTJ It 
is seen that for increasing R the deviation of the PFA value from the exact one decreases much faster 
than that for the energy. Already for R = 6, £ = 0.1, it amounts only a few percent. This allows for 
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TABLE II. The temperature dependent part of the force from the exact formula and from the PFA for e = 0.1 
and T — 1. 



R 


0.5 


1 


6 


h 


0.14 


0.56 


7.9 


/•PFA 
JT 


0.047 


0.31 


8.2 



the conclusion that at low temperature the PFA does not give correct results. However, for increasing 
temperature, the deviation of PFA from the exact results gets smaller and already at the border between 
low and medium temperatures the PFA gives quite good results for the temperature dependent part of 
the force. 




FIG. 2. The temperature dependent part Tt of the free energy (in arbitrary units) of a scalar held for e = 
and T = 1 as function of R calculated from the exact formula (solid line) and the corresponding proximity force 
approximation (dashed line). 




FIG. 3. The temperature dependent part Tt of the free energy (in arbitrary units) of a scalar held for e = 0.1 
and T = 1 as function of R calculated from the exact formula (solid line) and the corresponding proximity force 
approximation (dashed line). 

Now we consider the region of medium and high temperature by keeping d and T fixed while e — > 0. We 
start from representation ((T|) of the free energy for the electromagnetic field. We rewrite the Matsubara 
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sum by inserting unity in the form 

1= f° (55) 

J — OO 

Interchanging the orders of summation and integration we come to the representation 

rp pOO OO 

J7 =J * E ^-Cn)Trln(l-M(0). (56) 

Now we use the relation 

oo 1 OO 



2ttT 

71= — OO 



which is equivalent to the Poisson resummation formula, and the free energy takes the form 

T = \ r E ei</T Tr ^ (1 - M(0) . (58) 



2 7-oo 2tt 

^ n— — oo 



This expression differs from the vacuum energy ^ only by the sum and the additional exponential 
factor. For this reason we can apply the same methods as in [(| Q to calculate the asymptotic behavior 
for e — > 0. First we expand the logarithm and come to a representation in parallel to (fi2"1) . 



'-a/.i E »" £/WT> E^r E E II E W| ). <»> 

v n=-oo s=o m=-ooi=| m | V j=l n J= |m|-i/ 



where Z is given by Eq . (|43| . In (|59|) . in addition we substituted £ — > After that, Z(£/R) is a 

function of e only. Now we consider the limit e — > 0. For decreasing e, the main contribution to the 
integral over £ and and to the sums over the orbital momenta come from higher and higher values. 
Therefore we substitute the orbital momentum sums by corresponding integrals and change the variables 
according to 

t i r t tr 

t = - yl-T , 1 = -t, m=\ — u, 

e e V e 

[it 

n l = \ —n i (i = l,...,s). (60) 
V £ 

In this way we get an asymptotic expansion in the form 



T=— -V— I" dtte- 2t ^ [ dTT V 



= intVl-T a /(£iir) 

s=0 ' " " u *u v •■■ ■ n=-oo 



d M e-^ 2 (^) \ \\p^\2,-. (61) 



where 



Z&S - n (ff M i + na + n i+1 J e ^)+W( S+ i) , (62) 

with 771 = $3|_ (fLj — ™id-i) 2 - We stress again that these formulas are in complete analogy to the 
corresponding ones in [H, Hj] . The whole difference is in the sum over n and the T-dependent exponential. 
For instance, the factor Z as is the same as without temperature. Therefore we know it has an expansion 

Z as = l + 0(e). (63) 
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The corrections 0(e) depend only on the geometry and the kind of fields, but not on the temperature. 
Now we are interested in the pure temperature dependent contributions. Hence we drop these corrections 
and take Z &s = 1. After that J 7 , Eq. ff6T|) ■ can be simplified. First of all, one can perform the integrations 
over hi and /i, 



73 °° 1 />CO pi 1 

Now the integration over t is trivial. Also the integration over r can be carried out leaving the double 
sum 

R °° 1 °° e 2 T 2 

T= ~I^CP ^ (s+l) 2 ^ n 2 + 4 £ 2( s + 1 )2 T 2- ( 65 ) 



n— — CO 



Here the contribution from n = is just the zero temperature contribution, i.e., the vacuum energy. For 
the remaining sum over n^O we use the formula 



and the free energy can be written in the form 



V 1 = Q7rcotn ( a7r ) - 1 

^ a 2 + n 2 ~ 2a 2 [ ' 



T=-^{l + h(2dT)), (67) 



where we defined 

coth(m7ra;) 



h(x) = 90x 4 



(mnx) 3, (mirx) 4 



In the last step we substituted s + 1 — > m. 

In (|67p . the function h(2dT) describes the temperature contribution to the free energy for medium and 
high temperature at small separation. It was derived from the exact formula ([IJ . Now it is interesting to 
observe that this function, h(x) as given by Eq. (|6"5)) . is just the same as the function h(x), Eq. (P5)) . found 
in PFA in the preceding section. This can be shown, for example, by inserting the sum representation 
for g(x) into (|55|) . For this it is useful to rewrite (fl())) in the form 

g(x) = 45x 4 -|- r_£2__(_l + . (69) 

^yoyx y y J\y= m ^x 

Obviously, the integration in can be carried out explicitly and one comes just to Eq. (|6"5|) . In this 
way, we confirm the PFA (|36[) and with it all subsequent discussions in section 3. 

A similar statement holds for the force. In taking the derivative with respect to the separation d of 
(f6"7} . we have to consider the derivative of h(x)/x 2 . For this derivative the following formula holds, 

d h(x) 2 , 

~-Z3[9(x), (70) 



dx x x 3 

where g(x) is the same function as in the temperature dependent part of the free energy for parallel 
plates (fH| . This can be seen directly by inserting ([B"5f into the left side of ([70]) and comparing with (jTHJ) • 
In this way, also Eq.(pO| for the force in PFA is shown to coincide with the exact expression obtained 
from the T-matrix representation (JXJ) . 

Finally we mention that the function h[x) 1 like g(x), Eq. (ll7l) . has an inversion symmetry, 

i)=4-^, (") 

XI x z x^ 



which can be derived from (fTT)) and ([70]) . 
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V. THE FREE ENERGY AT LOW TEMPERATURE 



In this section we consider the low temperature expansion of the free energy, i.e., we assume T -C 1/R 
and T <C 1/d not restricting the relation between d and R. We start from representation (|4"4")l . The 
lowest approximation for T —¥ is, of course, the vacuum energy. The temperature corrections come 
from expanding Ft, (|45|) . for small T. In this section it is useful not to expand the logarithm and we 
take the representation 

Tt = ^ J di n T (0 i Tr [ln(l - M(*0) - ln(l - M(-i£))] ■ (72) 

For T — > 0, because of the Boltzmann factor, we can expand the remaining part of the integrand in ([721 
in powers of £ . Each additional power in £ will add a power in T. So we need the lowest odd power of £. 
The lowest even power of £ does not contribute to the difference between the logarithms which in fact 
represent the jump across the cut the logarithm has. 

First we expand M;^/(^), ((3|). We use the ascending series of the modified Bessel functions (see, for 
example [Hj]), 

'^=(S)"wtt y ^-^'U^if^-i^mU)- (73) 

where i v has only even powers of z and we continue to use the convention v = I + 1/2. Inserting these 
expansions into ([3]) and remembering that only I" give non-zero contributions for which I + V — I" is even, 
we observe that the first odd power of £ comes from the second term in K v (z) in (f75|) . Moreover, the 
lowest odd power comes from the lowest orbital momenta I and V only. For this reason from the orbital 
momentum sums only one term remains and the logarithm can be calculated directly. For Dirichlet resp. 
Neumann boundary conditions on the sphere we obtain in this way for the scalar field 



The logarithms are 



<o(0=-^ 2 + Ji? 3 e 3 + .... (74) 



ln(l-M D o(0)=ln(l-#)+i?e 



2L 

ln(l-M N o(C)) = ^e 2 -^ 3 + ... (75) 



and the jumps become 



: ,Tr [ln(l - M D (i£)) - ln(l - M D (-^))] = -2B£ 



iTr [ln(l - M N (z£)) - ln(l - M N (-z£))] = -|-R 3 £ 3 + . . . . (76) 

Finally we have to insert this into (1721) . The remaining integration can be done in terms of zeta functions 
and we come to 

J -D = _C(2) iiT 2 + ___ ) 
7T 

J -N = _2C(4) E 3 T 4 + _ _ _ (?7) 
7T 

These formulas are for Dirichlet boundary conditions on the plane. For Neumann boundary conditions 
on the plane we have to reverse the sign in the logarithms in (f?2"j) . A simple calculation in parallel to the 
above one results in 

n C(2) 2L-R 9 
j-n = ^ 7^4 + _ (78) 
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We see that the leading order in the low temperature expansion follows the behavior for low momentum 
of the corresponding T-matrix. For Dirichlet boundary conditions, from the s-wave, the lowest order is 
£ resulting in correction of order T 2 . For Neumann boundary conditions, due to the derivatives, this 
contribution is absent and the expansion starts with T 4 . 

For the electromagnetic field we have to take care of the polarizations. However, the off-diagonal 
entries of the matrix M start with an additional power in £ resulting from A. Since these enter the trace 
in quadratic combinations only, an additional power of £ 2 as compared to the diagonal contributions 
results. As a consequence, the off-diagonal entries do not contribute to the leading order at small T. 
So we are left with, separately, the contributions from the TE and TM modes. For these, the same 
considerations as for the scalar field apply with the only difference of (fit))) in place of ((!]) and the orbital 
momentum sum starting from the p-wave. We write down explicitly the relevant contributions to Mi y (£) 
for both modes and for relevant values of the azimuthal index m, 

wTEm=l/« & R 3 L 3 R2 \,2 R 3 f, R 3 \ ^ 

M "' ®=1^ + Sl( 1+ 10I?)* -t( 1+ 16PJ^ ' + -' 

^(0 = ^-1(1^)^^(1^)^..., 

<^)^ + S(-^)^?( 1 -^)^.... 

From here we calculate the jumps just like in (1761) . Introducing the symbolic notation A for these we 
note 

A TE > m=0 = --i? 3 £ 3 + ... , 
3 

A TE,m=l 2 161/ 3 + R 3 33 

A = 316L3-i?3 i?? +■•■' 

A TM.m=0 4 AL 3 + R 3 3 3 

A™' m=1 = \ R 3 e + .... (80) 
These expressions we have to insert into Tt and to integrate over £. In the end we get 

ir \ 3 16L 3 — R 3 3 4L 3 ~ R 3 J y J 

This is the low temperature correction for the electromagnetic field. We see that it gives an order T 4 
contribution like for the forces acting between parallel plates. 

Now we can consider (fgTj) for small separation. We can put d — directly and get 

J T ^i?¥ + .... (82) 

157T 

This expression is, as expected, different from the temperature corrections in PFA, Eq. (|3"Tj) . It gives the 
low temperature (this is T -C l/R) correction beyond PFA. 

It is also possible to consider (|5Tj) for large separation, or equivalently, for small R. One gets 



T T ="-^> 1-^ + ... )R 3 T\ (83) 



6C(4) (,_&_ 
it \ Ad 3 

The leading order coincides with the corresponding low temperature expansion of Eq. (4) in Q . 



VI. THE FREE ENERGY AT HIGH TEMPERATURE 



The high temperature expansion can best be analyzed starting from the original Matsubara sum (fTJ). 
The leading order for T — ¥ oo is given by the contribution with n = 0, i.e., by the lowest Matsubara 
frequency. 
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We separate the contribution from n = 0, 

T = TF + Fi(Td, TR) (84) 

and the contributions with w / are collected in F%. The leading contribution Fo is proportional to 
T and it is a function of the dimensionless ratio (|12[) . The function F\ depends on two dimensionless 
combinations. 

We consider the function Fq in more detail. It is given by the formula 

Fo(e) = iTrln(l-M(0)). (85) 
From (j3]) and using (|73j) we get for the scalar field with Dirichlet boundary conditions on the sphere 

7i,rD , m _(R\ V^rq + r + i/2) l+v 

Mi,A°) - \2L) 21^ + 1/2)1^ + 3/2)^ ■ ( 86 ) 

Here we took into account that from the sum over I" in ([3]) only the term with I" = I + I' did contribute. 
For Neumann boundary condition on the sphere we have to account for the derivatives. These result in 
a simple factor and we can write 

M»(0) = ^M£,(0) (87) 

For the electromagnetic field we have in addition the terms mixing the polarizations. However, since 
Alji is proportional to £ these do not contribute. So we are left with the additional factor and the 
difference in the derivatives. We come to 

M™(0) = M£j,(0)A{+/\ M™ (0 ) = £±!m£, (0)A|+/'. (88) 

These expressions can be inserted into ([85]) and for any finite e = d/R ^ the function Fq can be 
calculated numerically. 

We consider the limiting cases of e — > oo (large separation) and e — > (small separation). For large 
separation we simply have to expand the M; i //(0) into powers of \ je and see that only the lowest orbital 
momenta contribute. We obtain 



M D (0) = |- + ... ) M N (0) = -2^ 



M™,(0)=2^) 3 + ..., Af™ (0)=4 ^ 3 + .... (89) 

The corresponding functions i^o we get from (|85[) in this approximation by multiplication with (—1/2). 
For the electromagnetic field, i.e., adding the two contributions in the lower line in (I89|) . the result 
coincides with the lower line in Eq.(6) in [3J. 

In the opposite limit of small separation, i.e., for e — > we are faced with the problem that the 
convergence of the orbital momentum sums gets lost. This problem is essentially the same as at zero 
temperature in the same limit. Even more, it can be treated by the same methods, i.e., by calculating 
the asymptotic expansion of Fg(e) for e — > 0. The first step in this procedure is to expand the logarithm 
in (|85[) and to substitute the orbital momentum sums by corresponding integrations, 



F o& = -£ TXT / dl / dm II / dn i\ Z ( 9 °) 

S=0 S + L J0 J -l y j=1 Jl-m 

with 

s 

Z = l[M l+na+ni+1 (0)- (91) 
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We mention that this is, up to the sums substituted by integrals, the (n = 0)-contribution in (|42l) . Next 
we have to make an appropriate substitution of variables, 

/ = -, m = \/-Li, ni = \/ — ni. (92) 

£ V £ V £ 

This is the same substitution as used in §, Eq.(8) and in Q, Eq.(28), here however restricted to r = 1. 
Therefore we can use the asymptotic formulas derived there. For instance, in leading order we have 

M l+n , l+n ,(0) = (^) V2 e -2t-(«-n') 2 -M 2 (1 + o{e)). (93) 

This expression is the same for both boundary conditions and also for the modes of the electromagnetic 
field. For the latter there is in this approximation also no mixing of the polarizations, for details see Q. 
Making in (|90|) the substitution (|92|) we come for the electromagnetic field with (|93|) to 

- ~h E t^t f dte ~ 2t(s+1) r ^ ^ (s+1) I TT / ^= I e ' m + ■■■ ( 94 ) 

with r)i — X^iLol" - * — fii+i) 2 . Carrying out the integrations and the sum we obtain 

F (e) = -^- + ... (95) 

We mention that this formula coincides with the corresponding PFA result (|4"0)) although we did not 
assume d <C R- 

In general, the contribution of the zeroth Matsubara frequency is equivalent to a theory with a di- 
mension reduced by one. However, thereby it is usually assumed that the original theory (taken in its 
Euclidean version) has a symmetry between the spatial and the time dimensions. One must bear in mind 
that in the given case this symmetry is broken by the boundaries. 

Finally we consider F± in (|84|l . i.e., the contributions from the non-zero Matsubara frequencies. Here 
we observe for large T an exponential decrease of the functions Mu>(0), 

M M /(0) ~ e - 4 ™ Td (96) 

(d = L — R) such that F\ is exponentially suppressed for all d > 0. This is the same property as we 
observe in the case of parallel planes. This suppression holds for for ^ <C d, A <K R. If in addition 
d <C R, i.e. y <C d <C R, holds this suppression disappears. In that case, as before, higher and higher 
n and orbital momenta contribute. In fact, this is the limit RT — > oo in the region of high temperature 
as defined in section 3, Eq. pTj) . In that case it is appropriate use representation (|55|) and to proceed as 
it was done in the second part of section 4. As a result, F\ would deliver a correction beyond what is 
displayed in Eq. (|67|) . 



VII. CONCLUSIONS 



In the foregoing sections we calculated the free energy for a sphere in front of a plane at finite tem- 
perature for conductor boundary conditions. First we considered the case of small separation. Here we 
used both methods, the PFA and the functional determinant representation and considered both, the 
free energy and the force. In the region of low temperature, the thermal contributions to the free energy 
and to the force are very small and PFA does not hold for them. In the region of medium and high 
temperatures, which is the temperature region of experimental interest, we have reproduced the PFA 
from the exact method. For instance we have shown that the rule (f2"5"j) does hold in this case. 

We would like to mention that in [l3| exact finite temperature results were obtained for the interaction 
between a plane and a semitransparent curved surface described by a delta-potential. For weak coupling 
it was shown that the free energy of a scalar field coincides at all temperatures with the proximity force 
approximation corresponding to this geometry. 

Next we considered low temperature without restriction to small separation and derived the corre- 
sponding expansions for the free energy. Here the behavior of the thermal contribution to the free energy 
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depends on the boundary conditions and on the kind of the field. An interesting question on the in- 
terplay between temperature and geometry was raised in [3, [1[ . There it was observed that in an open 
geometry at low temperature lower powers of the temperature may appear as compared with the parallel 
planes. In general, we cannot support this conclusion from our calculations. Only for a scalar field with 
Dirichlet boundary conditions we observe Tt ~ T 2 , Eq. (|77p . For Neumann boundary conditions and for 
the electromagnetic field we see Ft ~ T 4 , Eqs. ([77|l (lower line), ([78)) and (|82|). 

In section 6 we considered the limit of high temperature. As expected, here the dimensional reduction 
works and the leading order contribution comes from the zeroth Matsubara frequency. At small separa- 
tion, the asymptotics could be calculated using the same methods as for zero temperature and the result 
coincides with that of PFA. 
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